A note on the QCD evolution of generalized form factors 
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Generalized form factors of hadrons are objects appearing in moments of the generalized parton 
distributions. Their leading-order DGLAP-ERBL QCD evolution is exceedingly simple and the 
solution is given in terms of matrix triangular structures of linear equations where the coefficients 
are the evolution ratios. We point out that this solution has a practical importance in analyses 
where the generalized form factors are basic objects, e.g., the lattice-gauge studies or models. It 
also displays general features of their evolution. 
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Generalized parton distributions (GPDs) (for nota- 
tions and discussion see the extensive reviews [l], 0, 0, 
0, H, [S, 7, 8, 9] and references therein) carry very rich in- 
formation on the internal structure of hadrons. In partic- 
ular, moments of the GPDs in the X variable, according 
to the polynomiality feature, can be written as polyno- 
mials in the ^ variable, assuming the form 

/ dX + l F^{X, i) = 2 ^ A2,+2.2^{t)e\ (1) 



with j = 0, 1, ... . Here X is the average fraction of the 
target's momentum carried by the struck quark, and ^ 
is the fraction of the momentum along the light cone 
passed to the target. We use the so-called symmetric 
notation, where X S [—1,1] and ^ G [0,1]. Indices ns 
and s denote the non-singlet and singlet distributions. 
The singlet GPD consists of the quark and gluon parts, 
F'{X,^,t) = (F^'«(X,e,i),XF^'^(X,e,t)). The gener- 
alized form factors Ank (GFFs) for the quark operators 
are the matrix elements between states i and / of the 
form (we take the even-parity operator and a spin tar- 
get) 

(,/(p')l?(0) 1^" iD^HD"' . . . iD"'^-'^ i^{0)\i{p)) = 
2p{^P^i...P^'-i>A,o(t)+ (2) 



k=2 
even 



with n = 1,2,3 k = 0,2, ... ,n, P = {p + p')/2, q = 
p' — p, t — q^, and 'ip denoting the quark field. The sym- 
bol D is the QCD covariant derivative, D = ^{D — D), 
and {. . . } denote the symmetrization of indices and the 
subtraction of traces for each pair of indices p-i . . . /in-i- 



The factor of 2 is conventional. A similar expression 
can be written for the singlet gluon GFFs [7[. 

For the simplest case of a spin-0 target with positive 
charge such as e.g. 7r+, Aio{t) is the charge form factor, 
while A2o(t) and —A22it) are the quark components of 
the gravitational form factors [l^, El ■ 

The GPDs undergo the QCD evolution with a change 
of the renormalization scale. Unlike the parton distri- 
bution functions (PDFs) or the distribution amplitudes 
(DAs), it is non-trivial to pass to the space of moments 
where the evolution is diagonal and then invert the trans- 
formation. This issue makes the case different and more 
complicated already at LO from the case of the PDFs, 
where the Mellin moments are used, or the DAs, where 
the Gegenbauer moments diagonalize the evolution. The- 
oretical tools have been developed to achieve the task, 
such as the Shuvaev transformations [l^, the dual rep- 
resentation of the GPDs [H, [13] or techniques based on 
the conformal moments cornplemented with the Mellin- 
Barnes transformations [l^ |3] ■ One r nay also solve the 
DGLAP-ERBL equations numerically [l3,|Il. Such ap- 
proaches are needed, if the whole GPD is demanded. 
Frequently, however, one is only interested or has ac- 
cess to a limited number of GFFs. The purpose of this 
note is to point out that in fact the LO DGLAP-ERBL 
evolution of the GFFs is inherently much simpler and 
straightforward to implement in practice, without any 
need of complicated mathematical transformations. The 
result discussed in this paper is implicitly present in nu- 
merous works concerning the evolution of GPDs (0] and 
references therein), but nevertheless we find it practical 
to present its explicit form, useful for those dealing with 
the GFFs only and not the full GPDs. 

Our starting point is the work of Kivel and Mankiewicz 
(T9I [20|. w hich elaborates the formalism of Balitsky and 
Braun [21[ on the QCD string operators in the coordinate 
space. We use Eqs. (18) and (20) from Ref. [3 for the 
LO DGLAP-ERBL evolution of the GPDs from the scale 
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fj,o, where one assumes they are known, to the scale fi: 



1 



IT \f3^ 



n=0 



+ 2n X 



UJ 



where 



dujF''%uj,^,t;fi'^) cos(w/3), 



dLuF'{uj,^,t;n^) sin(w/3). (4) 



For the non-singlet case the quantity L2n+i denotes the 
evolution ratio. For the singlet case £2n+2 forms a 2- 
dimensional matrix in the quark-gluon space. Explicitly, 



L2n+1 — 

r„ = 



^(^2) X 7-/(2/^0) 

7,?^ 7„^^ j ' 



(5) 



where the 7's denote the appropriate anomalous dimen- 
sions. Our further procedure is based on the observation 
that T{f3,£,,t; ^'^) is the generating function of the the 
GFF's. We then expand Eq. piH)) in /3 around using 
the series expansions 



i-iy 



1=0 



ll{m + iy. \2 



2:\ 2/+m 



(6) 
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forfc = 0,l,2,. 

Aio 
A32 

^54 



^30 
^52 



and Z = 0, 1, . . . , fc, or, explicitly, 
1 



^32 ^ 



(9Li - UL3 + 5L5)Al 
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+ ^{L3~L,)At,+L,A'i,, 



LsA" 



3^30' 



';{L3-L5)A%+L5Al^ 



^^52 1 



^50 = 



^5^501 



(9) 



where we have grouped the equations in the growing dif- 
ference of the indices n and i in Ani ■ The ellipses denote 
equations with n > 5. Since Li = 1, the vector form 
factor, of course, does not evolve. All other form factors 
in Eq. ^ change. While the standard form factors Ano 
retain their shape, i.e. Anoit)/Ano{t = 0) is not altered 
by the evolution, other genuine generalized form factors 
involve linear combinations and both their value at t = 
and their shape do change. Analogously, for the singlet 
case 



A2k+2,2i = r(2A: + 2) ^ (4m + 5)L 

m=0 

22J-2fc-i(_iyn-jp (j + m + I) AO 



k 

2m+2 

j=k-l 



(10) 



(i+l),2(i-fe+0 



r(2j + 2)r(m - i + \)T{k -m + l)T{k + m+l)'' 
for fc = 0, 1, 2, . . . and Z = 0, 1, . . . , + 1, or, explicitly, 

^22 = -^2^221 

3 

A44 — -(^2 — -^4)^22 + -^4^441 

Am = J^(11L2-18L4 + 7L6M^2 

+ i^(L4-L6K4 + L6<, 



where (.)„ denotes the Pochhammer symbol, and finally 
use the polynomiality property ([T]). As a result, polyno- 
mials in fi and ^ are obtained on both sides of Eq. (|3l3p . 
For subsequent values of the powers of (3 we compare the 
coefficients of powers of ^. As a result, the equations for 
the form factors from Eq. (|ll2p follow immediately. We 
use the short-hand notation Ank = Ank{X,S,,t] fi) and 
^nfe ~ Ank{X,$^,t; ^q). For the non-singlet case 



A20 

Aa2 
A&A 



^2^20, 

3 

-{L2 - -^^4)^20 + L4AI2, 
J^(llL2-18L4 + 7L6M^o 

+ Yy(-^4 ~ -^6)^42 + -^6^645 



^2^+1,2/ = kT{2k) (4m + 3)L2m+l Y 
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r(2j + l)r(m -j + l)T{k - m -M)r (fc -t- TO -I- I) 
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60 



-t^4^40 5 

— (-^4 - ^6)^40 + ^6^62! 



(11) 
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Note the identical structure of the first two groups in 
the above equation. Again, the shape of the form fac- 
tors Ano does not change. The sets of equations (|9|lip . 
although implicitly present in schemes involving the con- 
formal moments, have not, to our knowledge, been writ- 
ten explicitly and their practical importance has not been 
recognized. Since all quantities on the right-hand side are 
known, from any practical point of view the problem of 
the leading-order DGLAP-ERBL evolution of the GFFs 
is solved. 

Expressions for higher values of n may be obtained 
from the general expressions, however, it is the lowest 
form factors which are most relevant, as they can be 
obtained in the Euclidean lattice studies for the pion 
[H, H i, m and the nucleon [H, [11, [13, [111 . Equations 
(|9|lip are useful for the evolution of higher GFFs which 
eventually will be measured on the lattice as the accu- 
racy is increased, as well as for various model calcula- 
tions, where the results need to be evolved in order to 
compare to the data [29l |. 

Sequences of equations in (|9|lip separated by ellipses 
form (infinite) triangular matrix structures. Mixing oc- 
curs between the ns or s form factors where the difference 
n — z in Ani is fixed, for instance Aiq, A32, ^54, etc. This 
corresponds, according to Eq. ([2|), to the mixing of the 
t-channel states of the same angular momentum. Indeed, 
n — i is the number of Lorentz indices of the f-channel 
momentum q. These triangular matrix equations may be 
diagonalized, yielding the combinations 
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etc., which evolve autonomously with L„. The coeffi- 
cients in the above combinations are proportional to the 
coefficients of the Gegenbauer polynomial Cn^^(z). In 
fact, we are simply recovering the well known fact [X] 
that the conformal moments of the GPDs, 



r(3/2)r(n + l) 
2"+ir(n + 3/2) 



dXCCfJ'iX/OF^'^'iX, ^, t), (13) 



evolve autonomously at LO. The point is, however, that 
the diagonalization of Eq. p2)) is not necessary for the 
evolution of the GFFs, as for any practical purpose one 
can simply apply Eq. (I9|lip . 

Asymptotically, as /i^ — > 00 we have (for the positive- 
parity operator) Li = 1 and L„ ^ for n > 1. Hence, 
in the non-singlet channel 
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(14) 



In the singlet channel we give the sum of the quark 
and gluon components, whose second moment for the as- 
sumed vector case is related to the momentum sum rule. 



Denoting Snk = A'^^j^ + A^/. we find asymptotically 

'5'22 



D221 '-'44 y-'22i '-'66 ^'^22: ■ ■ 
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(15) 



Equations (|14I15|) comply to the following LO asymptotic 
forms of the quark and gluon GPDs, 



pns 
ps,q 



0{f - X"' 



1 - 



X^ 



Aw{t), 



e 

d{e-x^){S2,{t)+es22{t)) 

15 Nf X ( X^ 

e 



1 



XF 



s,G 



4^2 ACf + Nf ^ 
9{e-X'){S2o{t)+eS22{t)) 
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4Cf 



16^ 4Cf + Nf 



1 - 



±1 

e 



(16) 



The lowest form factors determine these expressions. 
Note, however, that the gravitational form factors S2o{t) 
and S22{t)) need not be equal. Only for the special case 
S2o{t) = ^S22{t) = S{t) one recovers the typically writ- 
ten form with the common factor (1 — ^'^)6{t). Since 
S22 corresponds to the coupling of a scalar, and 6*20 
to the traceless rank-2 tensor, there is no reason why 
S2o{t) = 5*22(0 should hold. This issue is related to the 
lack of factorization of the i-dependence in GPDs. 

In conclusion, we state that the generalization of the 
result to other channels is straightforward. For other 
probing operators one needs to simply use appropri- 
ate anomalous dimensions. For various targets (pion, 
nucleon) where different tensor couplings appear, one 
evolves the form factors separately for the independent 
structures. 

Gravitational and higher-order GFFs may be obtained 
from chira l q uark models of the GPDs for the pion 
[2l[33,[MI,[3l[3i,[3l[3i,|3i and for the nucleon 



A related quantity, the pi on-p hoton transition distribu- 
tion amplitude [33, '40', '41, '42*1, has also been obtained in 
quark models Refs. [43, 44, 45, 46, 47] and its moments 
undergo the QCD evolution in a similar way. Dynami- 
cal calculations test the simplifying but a priori unjusti- 
fied assumptions made in many phenomenological stud- 
ies. In particular, the widely assumed factorization of 
the ^-dependence is disproved. On the other hand, the 
reference scale, /^o, turns out to be very low in chiral 
quark models, around 320 MeV for the local models, and 
hence QCD evolution to experimentally accessible scales 
implies a long distance evolution. This is a case where 
the dilatation covariance may prove crucial, as it ensures 
the integrability of the renorm-group equations, hence 
the evolution path independence between two different 
scales [48*1. 

Actually, much of the explicit simplicity of the GFF 
evolution of Eq. (|9llip is linked to the LO approxima- 
tion and the conformal invariance of the evolution which 
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represents faithfully the dilatation group, a feature not 
automatically guaranteed in the current factorization ap- 
proaches at NLO, possibly inducing systematic errors. 
Renorm-group improvement of the GPDs might be im- 
plemented as previously done for the DGLAP evolution 
of the PDFs [43] as well as the ERBL evolution of the 
PDAs [4^ . While these NLO complications prevent writ- 
ing down a handy analytic solution for the GFF evolu- 
tion, the problem can be reduced to a set of coupled 
differential equations. It remains a tractable and much 
simpler alternative when a reduced set of GFFs is avail- 



able. See also Refs. 
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